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ABSTRACT 


We present results from a survey of a 1300 arcmin 2 region of the Orion B 
South molecular cloud, including NGC 2024, NGC 2023, and the Horsehead 
Nebula (B33), obtained using the Submillimetre Common-User Bolometer Array 
(SCUBA) on the James Clerk Maxwell Telescope. Submillimeter continuum ob¬ 
servations at 450 prn and 850 pm are discussed. Using an automated algorithm, 
57 discrete emission features (“clumps”) are identified in the 850pm map. The 
physical conditions within these clumps are investigated under the assumption 
that the objects are in quasi-hydrostatic equilibrium. The best fit dust tempera¬ 
ture for the clumps is found to be XJ = 18 ±4 K, with the exception of those asso¬ 
ciated with the few known far infrared sources residing in NGC 2024. The latter 
internally heated sources are found to be much warmer. In the region surrounding 
NGC 2023, the clump dust temperatures agree with clump gas temperatures de¬ 
termined from molecular line excitation measurements of the CO molecule. The 
bounding pressure on the clumps lies in the range log(/c _1 P cm 3 hU 1 ) = 6.1±0.3. 
The cumulative mass distribution is steep at the high mass end, as is the stellar 
Initial Mass Function. The distribution flattens significantly at lower masses, 
with a turn-over around 3 - 10 M 0 . 

Subject headings: infrared: ISM: continuum - ISM:clouds - ISM: individual (NGC 
2024, NGC 2023, B33, Orion B South) - ISM: structure - stars: formation 
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1. Introduction 


This paper continues an effort to quantify the necessary pre-conditions for star formation 
via observations of the continuum emission from cold dust at submillimeter wavelengths. 
Snch observations allow the determination of the mass spectrum and physical characteristics 
of the population of cold (typically 10-30 K), dense concentrations of dust, or “clumps”, 
each of which is expected to be a site of active, or eventual, star formation. It is also well 
known that stars form in groups (see e.g. Clarke et al. 2000, Elmegreen et al. 2000) and 
that these in turn originate from a hierarchy of dust clumps for which the mass spectrum is 
closely similar to the stellar initial mass function (IMF; see e.g. Motte et ah 1998; Johnstone 
et ah 2000b, hereafter Paper II; Johnstone et ah 2001, hereafter Paper III; Motte et ah 
2001). The recently launched Spitzer Space Telescope will observe many star-forming regions 
in the infrared and thus provide information on which dust and gas clumps contain deeply 
embedded protostars. Together, the submillimeter continuum dust maps and the infrared 
images will yield statistics on the number of prestellar, Class 0, and Class I protostars and 
significantly enhance our understanding of the lifetime in each stage and the changes in natal 
environment with protostcllar phase. 

Previous papers in this series on large area mapping at 850 /jrn have discussed image 
reconstruction techniques (Johnstone et ah 2000a, hereafter Paper I), and results for the 
Ophiuchus cloud (Paper II) and the northern part of the Orion B region (Paper III). In the 
present paper we turn our attention to the southern part of the Orion B molecular cloud 
(Lynds 1630). This region includes a number of objects well known to both amateur and 
professional astronomers. In particular, the Horsehead Nebula, B33 (see, e.g., Pound et ah 
2003), is seen in silhouette against the bright background of the HII region IC434 near the 
three bright stars of Orion’s Belt. At a distance of approximately 400 pc [Anthony-Twarog 
(1982) measured a distance of 390 pc] the Orion B South region also contains the young 
star clusters NGC2023 and NGC2024, and as part of the greater Orion Molecular Cloud is 
well-known as one of the nearest extensive regions of active star formation. An optical image 
showing most of the region discussed in this paper, and its relationship with the ionizing 
stellar system a Orionis, is given by Abergel et ah (2002). 

One of the first systematic searches for dense molecular cores in the Orion B region 
was carried out by Lada et ah (1991a; hereafter LBS) using the CS J = 2-1 transition, for 
which the effective critical density is about 10 4 cm 3 (Evans 1999). With a bearnwidth of 
1.8' a total of 42 dense cores were found to be concentrated in two relatively small regions 
of the Orion B molecular cloud; 19 of these lie within the region surveyed in the present 
paper. Seven of the latter objects were subsequently mapped by Launhardt et ah (1996) 
at 1300 fin l wavelength with an angular resolution of 12", comparable to that of the present 
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observations with the JCMT. From these surveys it is evident that active star formation is 
confined to two rather limited areas within the Orion B region. 

The region has been surveyed at many wavelengths. Early 2.2 /mi observations of a 
large part of the molecular cloud made by Lada et al. (1991b) identified four embedded 
star clusters, and subsequent widefield imaging of the stellar component is available via the 
2MASS images (Carpenter 2000). The distribution of relatively warm dust was mapped by 
IRAS, although for our purposes with limited angular resolution, and in far infrared emission 
at 138 and 205 /mi (Mookerjea et al. 2000). The molecular component has been surveyed 
by a number of authors, e.g. Miesh and Bally (1994) in the 13 CO J = 1-0 transition, and 
by Kramer et al. (1996) with lower angular resolution in the 13 CO and 12 CO J = 2-1 and 
3-2 transitions. 

The Horsehead Nebula in particular has been extensively investigated. Pound et al. 
(2003) have obtained an interferometric image in the CO J = 1-0 transition of the Horsehead 
Nebula with an angular resolution of 10". Abergcl et al. (2003) combine mid-infrared 
ISOCAM data in the LW2 and LW3 bands with CO, 13 CO and C 18 0 1-0 and 2-1 observations 
to model the transient heating of the outer edge of the dust cloud by photodissociating 
radiation from the relatively nearby 09.5V system a Orionis. Teyssier et al. (2004) have 
obtained detailed observations of simple cyclic and linear hydrocarbons, which they show to 
be distributed similarly to CO within the western edge of the nebula. Finally, early results 
from submillimeter-wavelength mapping of the Horsehead Nebula region have been reported 
by Sandell et al. (2001); these latter archival data have been included in the dataset for the 
present work, which covers a total area about 1300 arcmin 2 in extent. 


2. Observations and Data Reduction 

The data covering Orion B South were obtained between 1997 and 2000 using SCLIBA 
at the JCMT. The JCMT archive, hosted by the CADC 1 , was used to find all relevant 
observations toward the region. A total of 135 individual scan and jiggle maps were retrieved, 
which, when combined, cover a total observed area of 1300 arcmin 2 . Each individual 3-arcsec 
cell in the final data array was measured approximately 40 times at 850 /mi and 90 times 
at 450 /mi. Thus, the total time spent observing each cell was approximately 5 and 11 
seconds respectively. The weather conditions varied significantly during the many nights of 
observing. The mean optical depths at 850 /mi and 450 /mi were on average quite favorable, 


1 The Canadian Astronomy Data Centre (CADC) is operated by the Dominion Astrophysical Observatory 
for the National Research Council of Canada’s Herzberg Institute of Astrophysics. 
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7~85o = 0.26 ± 0.07 and T 450 = 1.4 ± 0.4. 

The raw bolometer data were reduced using the standard SCUBA software (Holland et 
al. 1999) to flat-field, extinction correct, and flux calibrate. The extinction corrections and 
flux calibrations as a function of time were obtained from the JCMT calibrations archive 
(Jenness et al. 2002). The extinction corrections are tabulated as a function of time based 
on a comparison of skydip measurements taken with the JCMT during normal observing 
and zenith optical depth measures taken every few minutes at 225 GHz by the CSO tau 
meter. The flux calibration values are tabulated from all calibration measurements taken 
during extended periods of time, e.g. semesters, over which the telescope and instrument 
were not known to have changed significantly. Typical uncertainties in the flux calibration 
can be determined from Jenness et al. 2002 (see their Figures 1 and 2 ). At 850 jum the 
uncertainty is approximately ±10 % while at 450 /mi the uncertainty is closer to ±25 %. 
These values are consistent with the uncertainties determined by individual observers using 
calibrator sources. Most of this uncertainty is due to changes in the telescope surface due to 
temperature and gravity during the night and resultant changes in the beam profile, which 
can be quite significant at 450 /mi. 

The data were further reduced and transformed into the final maps using the matrix 
inversion technique described in Paper I. This method has many advantages over the standard 
facility data reconstruction procedure, in particular the ability to use all data taken at 
the telescope, regardless of chop configuration or mapping procedure, and the ability to 
properly weight each individual measurement based on bolometer and sky noise. Significant 
low-amplitude, large-scale features remain after reconstruction and may be real; however, 
artificial features produced by weather variations and bolometer drift over time will also 
produce this type of structure since reconstruction of chopped data amplifies the longer 
spatial wavelengths (chop throws range from 20" to 65"). We filter out these large-scale 
features by convolving the map with a large Gaussian profile (a = 135") and subtracting the 
resultant smoothed map from the reconstruction. In order to minimize the effect of over¬ 
subtraction in regions with bright sources, all pixels greater than 0.25 Jybeam -1 at 850 /mi 
and 5 Jybeam -1 at 450 jum have been removed in the construction of the smoothed map. 
The final, reconstructed 850 /mi and 450 /mi maps, after filtering, are shown in Figures 1 and 
2 along with maps of the estimated uncertainty at each position. Note that the uncertainty 
in the measurement drops significantly in regions where multiple observations were made, 
such as the Horsehead region. 
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3. Data Analysis: Clump Properties and Distribution 
3.1. Identification of Clumps 

We have noted that size scales significantly larger than the maximum chop throw are 
unlikely to be correctly registered in the final images, and such structures have been removed. 
The data are thus amenable to an automated clump-finding analysis for structures of size 
scales up to approximately 130". 

Clumps within the Orion B South region were identified in substantially the same man¬ 
ner used for Ophiuchus (Paper II) and Orion B North (Paper III), allowing for a direct 
comparison of the results. The flattened 850 /jrn map was searched for structure using a two- 
dimensional version of the clump-finding algorithm, clfind (Williams et al. 1994). This 
algorithm divides the map into individual clumps by utilizing minimum flux boundaries, 
assigning each cell in the map above a limiting threshold to a particular clump. Thus, unlike 
Gaussian clump-finding algorithms, the regularity of a clump is not preset and clumps are 
allowed to have arbitrary shapes. A comparison of this clump-finding algorithm and “by 
eye” techniques is discussed in Paper III, where it is shown that the two techniques produce 
very similar results; however, with no observer bias in the present case. 

The typical sensitivity in the raw map is similar to that found in Paper III, about 
40 mJy per pixel, or 10 mJy per beam at 850 //m. Comparison can also be made with the 
millimeter-wavelength maps of Orion B produced by Launhardt et al. (1996) which had a 
sensitivity of approximately 130 mJy per 30" beam. Assuming the dust opacity scales as 
v 2 the mass sensitivity of the Launhardt observations is about sixteen times less good than 
in the present work. In order to maintain a uniform clump sensitivity across the map, the 
typical noise value in the flattened map, approximately 30 mJy per pixel, was used as input 
to clfind. 

Initially, 69 clumps were obtained by the clfind routine; however, twelve were removed 
due to location within the map, leaving 57 clumps. All removed clumps were near the edges 
of the map where residual noise is highest. Following Papers II and III, the size of each clump 
was computed by taking the area of each clump and determining the “effective radius” of 
a circle with equivalent extent. It should be noted, however, that a number of the clumps 
deviate significantly from circular symmetry. The properties of the 57 clumps are presented 
in Table 1 and the locations and relative sizes of the clumps are shown in Figure 3. The 
total flux for each clump has not been corrected for the effect of the telescope error beam, 
in keeping with previous papers in this series. The error beam at 850 /ini is at most ~ 10 %, 
much less significant than the calibration uncertainty. 
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An initial estimate of the mass for each clump was computed from the total flux within 
the clump boundary under the assumption that the measured flux was due to thermal emis¬ 
sion from optically thin dust particles (c.f. equation 1 in Paper III): 


A/clump = 0.59 X Si 
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/ 17K 
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^850 


0 . 02 cm 2 g _1 
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( 1 ) 


Following the results of van der Tak et ah (1999) we take the opacity per unit mass 
column density at 850 /mi to be ^sso = 0.02cm 2 g _1 , a factor of two higher than that used 
in our previous papers. The distance to Orion B South is assumed to be d $50 = 400 pc, 
somewhat less than the 450 pc distance used in Paper III for Orion B North. Anticipating 
the results of the detailed analysis below, the typical dust temperature for each clump is 
taken to be T d = 20 K, producing clump masses in the range 0.37 to 90 M 0 . The four most 
massive clumps are atypical however; they are coincident with known far infrared sources 
(Mezger et al. 1988) and have much higher dust and gas temperatures, T d > 50 K (Mangum, 
Wootten, & Barsony 1999; Fissel et ah 2005). For these four sources the derived masses are 
substantially decreased if the higher temperature is adopted, lowering the maximum clump 
mass found in Orion B South to less than 30 M 0 . The derived mass versus size relation for 
the clumps is plotted in Figure 7, assuming T d = 20 K for all clumps. 


In Paper III, the clumps in Orion B North were assumed to have a lower value of k 
than used here, and a higher typical dust temperature, T d = 30 K. As well, the distance was 
taken to be 450 pc rather than 400 pc. To compare derived masses between Paper III and 
this paper on a uniform scale, multiply by 0.7 the masses quoted in Paper III. 


3.2. Clump Distribution 

Within the Orion B South region the three objects B33 (also known as the Horsehead 
Nebula), NGC2024, and NGC2023 contain the large majority of the star-forming activity 
and submillimeter-wavelength emission clumps discussed in this paper. We present expanded 
images of these regions at 850 and 450 /jrn from the present work, in each case together with 
a corresponding optical image from the Digitized Sky Survey (DSS) and the 8 -/mi image 
obtained by the MSX satellite (Price et al. 2001). 


3.2.1. The Horsehead Nebula 

Of the three regions, only B33 presents aspects common to all wavelengths (Figure 4). 
The top of the Horsehead “neck” and “head”, which provide a sharp boundary to the dark 
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cloud silhouetted against the bright background of the optical emission nebula IC434, is 
revealed as a region of cold dust at submillimeter wavelengths. The western edge, heated by 
a Orionis further to the west (Abergel et al. 2003), appears bright at 8 /mi, and contains 
a compact feature at rnid-IR wavelengths which could a prestellar object. The latter is not 
differentiated at sub-mm wavelengths. 

The clump-finding routine identifies two 850 /mi sources (054089-02271 and 054091- 
02278; see Table 1), both parts of the “horse’s mane”, in this region, but neither are specif¬ 
ically associated with the compact rnid-IR feature noted above. Both clumps have moder¬ 
ately cool submillimeter-wavelength dust temperatures (Table 2). A third bright dust clump 
(054112-02273) 3' to the east, within the “horse’s neck”, has no obvious counterparts at rnid- 
IR or optical wavelengths, but is amongst the warmer sources based on the submillimeter 
data (see Table 2). Other submillimeter sources further to the east are parts of a clumpy 
ridge of emission extending south from NGC2023. Mid-IR emission possibly associated with 
this ridge is displaced slightly to the west, and likely results from heating by a Orionis. 


3.2.2. NGC2023 

Optical and rnid-IR images of NGC2023 bear no resemblance to the submillimeter- 
wavelength images (see Figure 5), and very little to each other. The optical image is dom¬ 
inated by a bright reflection nebula with some patchy absorption by foreground dust. The 
illuminating star for the reflection nebula is HD 37903; heated dust gives rise to the struc¬ 
ture shown the 8-/im MSX image and also in ISOCAM data (Abergel et al. 2002). The 
submillimeter-wavelength data however reveal a strongly clumped structure containing 12 
individual components. Comparison with the LBS dataset shows that three CS clumps lie 
within NGC2023. LBS 34 can be identified with the submillimeter clump 054162-02172, 
while 054142-02160 is likely to be related to LBS 36. Interestingly, no LBS clump is identi¬ 
fied with the brightest continuum source in NGC2023, 054141-02180; this object is known 
to harbor a very cold Class 0 source with a well-collimated molecular outflow (see Reipurth 
et al. 2004, and references therein). The remaining LBS clump within the field in Figure 5, 
LBS 39, is not clearly associated with compact submillimeter continuum emission. 


3.2.3. NGC2024 


As a major star formation region NGC2024 has been extensively studied. As shown in 
Figure 6 (upper right panel) the optical appearance is of an extensive ionized region overlaid 



with a dense foreground dust lane; the visual extinction toward the deeply-embedded object 
IRS 2 (Grasdalen, 1974) has recently been determined as > 27 mag at 2.2 /jrn (Lenorzer et al. 
2004). Although IRS 2 contributes significantly to the ionization, it was known for some time 
that this star did not meet the 09-09.5 stellar equivalent required. Recent work by Bik et 
al. (2003) suggests that a second fainter star, IRS 2b, only 5" from IRS 2, actually provides 
the necessary ionization. Centimeter-wave aperture synthesis continuum observations of this 
region (Barnes et al. 1989) show an HII region with a sharply-delineated “bay” to the south. 
This structure cannot be seen at infrared wavelengths, but is replicated in the mid-infrared 
emission at 8 /irn (lower left panel, Figure 6). 

The submillimeter wavelength observations of NGC2024 (see Figure 6) show a struc¬ 
ture consisting of two major clumps of dust emission bordered by a number of filamentary 
structures. Earlier millimeter-wave observations by Launhardt et al. (1996) reveal that the 
two major clumps of emission, as well as the CS clump LBS 33, coincide with this structure. 
At 633 M 0 , this is the most massive molecular clump found by Lada et al. (1991) in the 
Orion B South complex. 

The SCLIBA observations in this paper identify a total of 16 clumps within the region 
shown in Figure 6, including the four brightest sources in the map. These objects collectively 
and individually represent by far the most massive concentration of clumped dust within 
the entire Orion B South region. Notably, all of these objects have been assigned dust 
temperatures of 50 K, much higher than the mean dust temperature in the larger region, 
on the grounds that the 850 and 450 /irn fluxes are inconsistent with a modified blackbody 
model. Molecular line observations (Fissel et al. 2005) confirm that the clumps are warm. 


4. Modeling Clumps as Bonnor-Ebert Spheres 

The nature of the submillimeter-wavelength clumps found in nearby molecular clouds is 
at present quite controversial. Whether the clumps are long-lived, slowly evolving entities, 
perhaps the result of a gradual release of magnetic support (Mestel & Spitzer 1956; Shu 1983; 
Nakano 1984; Shu et al. 1987; Mouschovias & Ciolek 1999; Basil & Ciolek 2004), or transient 
turbulent structures (Scalo 1985; Mac Low & Klessen 2004 and references therein) has yet to 
be definitively settled. The recently launched Spitzer Space Telescope will play an important 
role by determining which clumps have embedded protostars. Regardless, molecular line 
observations show that the line widths of such dense cores are narrow, indicating that non- 
thermal motions do not overwhelm the dynamics (Barranco & Goodman 1998; Goodman 
et al. 1998), and implying that equilibrium models should reasonably reproduce the clump 
conditions. 
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In Papers II and III, physical conditions for the clumps were estimated through modeling 
the observations as Bonnor-Ebert spheres (Ebert 1955; Bonnor 1956; Hartmann 1998). In 
particular, the observable size R" e s, the flux S 850 , and degree of concentration C for a Bonnor- 
Ebert sphere map directly to the physical mass, temperature, and bounding pressure of the 
object. The concentration is found from the observable size, total flux, peak flux Esso, and 
the FWHM beamsize B" (c.f. equation 2 in Paper III): 


1.13 B" 2 S 8 5o 
TT R " e fp F 8 50 


( 2 ) 


A number of important assumptions are required when applying this technique to avail¬ 
able submillimeter-wavelength data. First, the conversion from flux to temperature requires 
knowledge of the dust opacity (see above) and the assumption that the dust and gas tempera¬ 
tures are related. Given the high densities for these clumps (n > 10 4 cm -3 ) it is reasonable to 
assume the dust and gas temperatures are similar and thus we take Tj = T g . In a departure 
from our earlier models, we assume also that an additional turbulent component of pressure 
exists within the clumps, equal to the thermal pressure at every location. This change in 
our analysis is motivated by observations by Goodman et al. (1998), where no clumps were 
found to have less than approximately fifty percent turbulent support using dense gas tracer 
molecules 2 . Given that the bulk of the molecular cloud has strongly non-thermal support, 
we take the transition case, equal contributions from both thermal and non-thermal, as an 
approximation for the clump support. Finally, determination of both the clump boundary 
(i.e. size) and degree of concentration (i.e. the importance of self-gravity) are complicated 
by the limited angular resolution and low sensitivity of the observations. Most of the clumps, 
however, lie well below the concentration at which clumps become unstable (see Figure 8) 
and thus are well represented by almost constant density models. For these cases, determi¬ 
nation of the physical properties, mass, temperature, and bounding pressure are not strongly 
dependent on the determination of the exact value of the central concentration. The clumps 
with the more extreme concentrations, located in the unstable part of the diagram, are 
marked with diamonds in Figure 8. The known far infrared sources lie among this group. 

Applying the Bonnor-Ebert sphere analysis used in Papers II and III, the physical 
conditions for each clump are plotted in Figure 9. Excluding the extreme concentration 
sources (marked by diamonds), the typical clump temperature is T c i = 18±4K. The pressure 


2 We note, however, that the level of turbulent support may be strongly dependent on the region. Recent 
results by Tafalla (private communication; see also Tafalla et al. 2004) show that in Taurus the fraction of 
pressure support due to turbulent motions is much less than fifty percent. 
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required to bound the clumps is log(/c _1 P cm 3 K _1 ) = 6.1 ±0.3. The scatter in T d and k~ l P 
values is only marginally larger than the uncertainty in any individual source measurement. 
A comparison with the results of our previous papers must take into consideration the change 
in the opacity law and the inclusion of additional turbulent support within each clump. 
Together these changes almost exactly cancel in the determination of the bounding pressure 
and approximately halve the computed internal temperature. With this in mind, the clump 
temperatures in Orion B South appear similar to those in Orion B North (Paper 111) and 
somewhat warmer than those in Ophiuchus (Paper 11). The bounding pressure is lowest in 
Orion B and highest in Ophiuchus. As noted in Paper III, this is in general agreement with 
the known conditions within the clouds; Ophiuchus has a significantly higher central column 
density, and hence internal pressure, than Orion, where observations have been obtained. 
In Orion the typical peak extinction through the cloud is 8 magnitudes with a maximum 
of 17 magnitudes (Lambardi and Alves 2001), whereas in the core of Ophiuchus the peak 
extinction reaches 36 magnitudes (Johnstone et al. 2004). Thus the Orion molecular clouds 
are bathed in an enhanced interstellar radiation held which warms the region. 


4.1. Clump Temperatures Derived from Submillimeter Colors 

In principle, the dust temperature within each clump can be estimated from the sub- 
millimeter spectral index, 7 , where S'(zj) oc z/ 7 . Assuming that the dust opacity follows a 
power-law with k { v ) oc A 5 , and that the dust emission can be approximated by a blackbody 
at temperature Td, the dust emission is given by 

S(u) oc A 3 B v {T d ). (3) 

In practice, however, determining the dust temperature is complicated by the fact that the 
dust opacity power-law coefficient, (5, is poorly constrained and often varies with location 
and physical conditions (Goldsmith, Bergin, & Lis 1997; Visser et al. 1998; Hogerhcijde & 
Sandell 2000; Beuther, Schilke, & Wyrowski 2004). As well, the 450 /mi and 850 /mi beam 
sizes at the JCMT are significantly different and the 450 //in beam is not well approximated 
by a single Gaussian function (Hogerheijde & Sandell 2000); at least half of the flux in 
the 450 /jrn beam is associated with a broad component about 30" in extent. Thus, simple 
convolution of the central 8.5" 450 //m beam with the 14.5" 850 /jrn beam will significantly 
overestimate the contribution of 450 /jrn emission on large scales. A sophisticated solution 
to determining f3 is to use a Fourier Transform Spectrometer to produce a low resolution 
spectrum of the 850 /jrn emission alone. At present, however, this is only possible for very 
bright sources (Friesen et al. 2005). Here we utilize the 450//m and 850/mi emission from 
SCUBA and follow the approach of Reid & Wilson (2005) in preparing the maps for analysis. 
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We convolve the 850 /jrn map with a model of the 450 /iin beam and the 450 /irn map with a 
model of the 850 /irn beam. For this purpose the 450 /irn beam is modeled as two Gaussian 
components with 8.5"and 30" FWHM, and relative peak intensities of 0.95 and 0.05. The 
850 /irn beam is modeled as two Gaussian components with 14.5"and 30" FWHM, and relative 
peak intensities of 0.95 and 0.05. 

Assuming that (3 = 2 (Hildebrand 1983), the central dust temperature for each clump 
was estimated using the total flux in the convolved 850 /irn and 450 /irn maps by inverting 

Ai = 24 H (ft) ~ 1 

^850 exp — 1 

The derived values are provided in Table 2 and a comparison with the derived Bonnor- 
Ebert model temperatures is shown in Figure 10. It is clear from the figure that while 
there is qualitative agreement on which clumps are warm, the spectral energy distribution 
derived temperatures for most clumps are poorly constrained by the available submillimeter- 
wavelength data. Some of the clumps, particularly at the faint end of the distribution, 
have 450 /irn to 850 /irn flux ratios which are inconsistent with modified blaekbodies having 
(3 = 2. In Table 2 and Figure 10 we assume that these sources are warm ( T,( = 50 K), 
although alternative explanations include (3 > 2 and/or measurement unreliability for the 
derived fluxes. We have already noted in Section 2 that the typical uncertainty in the flux 
determination at 850 /irn is 20% whereas at 450 /irn the typical uncertainty is 50% (Jenness 
et al. 2002), notwithstanding the complication of convolution discussed above. Setting aside 
the atypical high flux ratio sources, the mean dust temperature of the clumps derived from 
fitting the spectral energy distribution between 450 /irn and 850 /irn is Tj = 18 ± 6 K. 

Once protostcllar source catalogues obtained from Spitzer Space Telescope observations 
become available for the Orion B South region, a detailed investigation between the observed 
concentration and temperatures of clumps derived from the continuum observations and 
the presence of internal sources will be possible, yielding significant improvement in our 
understanding of the evolution of these prestcllar and protostellar envelopes. 

4.2. Clump Temperatures Derived from Molecular Line Observations 

Mapping of CO isotopomers was carried out at the JCMT during the same observing 
runs, when sky conditions did not permit continuum observing. Scan maps of part of the 
Orion B South region, NGC2023, were obtained in 13 CO J = 2-1, C 18 0 J = 2-1, and C 18 0 
J = 3-2. The maps are shown in Figures 11, 12, and 13. Figure 14 shows 850 /irn continuum 
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contours overlaid on the C 18 0 J = 2-1 map. The J = 2-1 maps were acquired using 7.5" 
sampling, which is about one third of the 21" beamwidth. The C 18 0 J = 3-2 map has a 
beamwidth of 14" and was acquired using 5" sampling. 

The 13 C0 to C 18 0 J = 2-1 line ratio provides the optical depth over the map, on 
the assumption of LTE. Goldsmith, Bergin, and Lis (1997) find that the abundance ratio, 
13 C0/C 18 0, lies in the range 6-15. In the present analysis we assume that the abundance 
ratio of 13 CO to C 18 0 is 10. As expected, the C 18 0 emission is optically thin, with r being 
typically less than 0.1, and ranging up to a maximum of 0.2. In addition to being optically 
thin, the C 18 0 J = 3-2 map has the same resolution as the 850 /jrn map and so is the most 
useful for a comparison of dust and gas emission. The C 18 0 and continuum maps agree 
quite well in general appearance, but a closer examination shows that the C 18 0 emission 
is somewhat smoother, less clumpy, and more extended than the dust continuum emission. 
These differences might be due in part to the removal of extended continuum emission by 
the SCUBA chopping technique, in part to different dependence of CO and dust emission 
on temperature, and in part to the freeze-out of CO molecules onto grains. 


We have obtained the CO excitation temperature from the ratio of the two rotational 
lines of C 18 0. The relevant equation is 


T P , 


16.59 log e 


4 7r(3 — 2) ' 
9/*(2-l)_ ‘ 


( 5 ) 


This equation is obtained from standard LTE and optically thin expressions and is derived 
in Mitchell (1993). Here, Ir is the integrated intensity of the appropriate spectral line, with 
the line temperature expressed in units of radiation temperature ( Tr = TX/r/mb , with the 
assumption that the beam filling factor is unity). We use appropriate beam efficiencies for 
the two transitions (r) m i, = 0.69 for the 2-1 line and 0.63 for the 3-2 line) and convolve the 
C 18 0 J = 3-2 observations to match the 21" beam of the J = 2-1 data. We End that, 
for most pixels over the map, the C 18 0 excitation temperature lies in the range 14 K to 
20 K, with a mean of 16 K. The (rare) extreme values are 10 K and 25 K. The uncertainty 
in individual values, due only to noise in the data, is ±2 K. These excitation temperatures 
are quite consistent with the average temperature of 18 ± 4 K found above for clumps in the 
NGC2023 region, using a Bonnor-Ebert analysis. 


Five of the continuum clumps in Table 1 fall inside the region of NGC2023 that we 
mapped in all three CO transitions. Using the continuum positions of these clumps, we have 
calculated the C 18 0 excitation temperature from the above equation and the C 18 0 column 
density using the expression given by Mitchell et ah (2001). Optically thin values of column 
density were multiplied by the factor r/(l — exp(—r)) to correct for optical depth. The 
C 18 0 column density was converted to a H 2 column density using the ratio N(H 2 )/N(C 18 0) 
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= 6 x 10 6 . In Table 3 we list, for these five continuum clumps, the temperature from 
the Bonnor-Ebert analysis, the dust temperature from the 850/mi/450/xm ratio, the C 18 0 
excitation temperature, the H 2 column density derived from the 850 /jrn peak flux per beam 
(assuming the BE derived temperature), the H 2 column density from the C 18 0 observations, 
and the ratio of the two column density determinations. 

The gas temperatures agree very well with the dust temperatures. This close agreement 
is rather unexpected in view of the large uncertainties in the two dust temperature techniques. 
The single exception is 054153-02189, where the dust temperature from the submillimeter 
continuum ratio is anomalously low. Column densities derived from CO and dust emission 
agree reasonably well for four of the clumps (054153-02189, 054162-02172, 054167-02171, 
054173-02170) but differ by a large factor for one (054141-02180). A similar analysis for 16 
clumps in Orion B North (Mitchell et al. 2001) resulted in masses from dust continuum 
exceeding masses from C 18 0 by factors of a few. The present closer agreement is due in part 
to the use of a larger grain opacity in this paper. 

Freeze-out of CO molecules on grains is a common phenomenon in cold cores (e.g. 
Jorgensen, Schoier, and van Dishoeck 2005) and could explain the high N(850)/N(C ls O) 
ratio for the one clump. However, since all five clumps have similar temperatures, why has 
molecular depletion acted in only one of the five clumps? Various uncertainties enter the 
column density determinations, including the value of the grain opacity and CO isotopic 
abundances, so that the ratio, N(850)/N(CO), could easily be increased by a factor of two 
or three and become consistent with some CO freeze-out in all cases. An examination of the 
C 18 0 J = 3-2 map shows that obvious clumps are present at only two of the five positions, 
namely continuum clumps 054141-02180 and 054167-02171. The apparent agreement in the 
case of the Erst clump is deceptive, however. The C 18 0 peak is displaced 20" to the southeast 
of the continuum peak. Although there is an extension of 850 /jrn emission from this clump 
towards the southeast, the clump-finding algorithm did not detect an independent clump at 
this position. Examination of the higher resolution 450 /jrn map does show a separate peak 
at the C 18 0 position. In fact, this second peak was detected by Sandell et al. (1999) in 
an early SCUBA observation. They called this clump 054141-02180 NGC 2023 mml and 
they named the source to the southeast NGC 2023 mm2. The position of the latter source 
coincides with the C 18 0 peak. From a variety of evidence, including a CO bipolar outflow, 
Sandell et al. (1999) suggested that NGC 2023 mml is a Class 0 YSO. It appears that strong 
CO depletion has occurred in mml but not in mm2. Clump 054167-02171 is bright in both 
continuum emission and in C 18 0 emission. It is almost certainly significant that this clump 
is the warmest of the Eve clumps, with a temperature of 24 K. This is warm enough that CO 
depletion is not expected to occur (e.g. Figure 1 of Jorgensen et al. 2005). 



4.3. The Clump Mass Distribution Function 


Stars form inside dense condensations within Giant Molecular Clouds. At least some of 
the clumps investigated above are expected to form, or to be forming, stars, and thus under¬ 
standing the physical properties of these clumps and their environment will help constrain 
formation mechanisms. One important bulk property of the clumps is the distribution of 
clumped mass. Following the pioneering study of the Ophiuchus molecular cloud complex 
by Motte et al. (1998), a number of submillimeter-wavelength surveys have been completed, 
each resulting in essentially the same conclusion, despite the objects found being in differ¬ 
ent environments. The mass function of the submillimeter condensations (Testi & Sargent 
1998; Paper II; Paper III; Motte et al. 2001) essentially mimics the stellar IMF where 
N(M ) oc M -1 ' 35 for stars more massive than M ~ 0.5 M 0 (Salpeter 1955) and the bulk of 
the mass resides at the low-mass end (Scalo 1986; Kroupa et al. 1993). This result is in 
direct contrast with the mass function of molecular clouds (Williams & McKee 1997; McKee 
& Williams 1997; Kramer et al. 1998; Williams, Blitz, & McKee 2000 and references therein) 
where N(M ) oc M~ a5 and the majority of the mass is found in the few large clouds. 

As in Papers II and III, the clump mass distribution is plotted in terms of cumulative 
number N(M ) of clumps with masses greater than M. Figure 15 presents the Orion B 
South data for both a constant flux-to-mass ratio (for TJ = 20 K) and for masses calculated 
using the Bonnor-Ebert sphere analysis. The horizontal dotted lines indicate the possible 
variation of the distribution due to uncertainties in the Bonnor-Ebert derived temperatures 
(as discussed in Paper III). The two dashed lines in Figure 15 denote power-law fits with 
N(M ) oc M~ a , where a = 0.5 (shallow line) or a = 1.5 (steep line) and are presented to 
guide the eye. The cumulative mass function determined using Bonnor-Ebert derived masses 
is steep, as is the stellar IMF at the high mass end, but flattens significantly around M ~ 3 
- 10 Mg. This is in contrast to the clump IMF measured in Orion B North (Paper III) and 
Ophiuchus (Paper II) using the same technique and with similar sensitivity. In both those 
distributions the steep mass function continued to M r^j 1 M q . 

One explanation for the observed discrepancy may be the environmental conditions in 
which the clumps are found. The typical clump in Orion B South is more extended, and 
thus less dense, than those found in Orion B North or Ophiuchus. The lower density of the 
clumps may imply a lower bounding pressure, in agreement with the Bonnor-Ebert analysis 
above. Why the clump mass distribution should be strongly affected by a relatively small 
change in bounding pressure is, however, unclear. The Jeans’ mass at which clumps become 
unstable varies only with the square root of the bounding pressure. More likely, the effect is 
due to the sensitivity limits of the survey. The mass versus radius distribution of clumps in 
Figure 7 reveals that incompleteness in clump identification may already become important 
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at M ~ 5 M 0 , explicitly because the clumps are more diffuse, lower in surface brightness, 
and thus harder to detect. 

There is significant discrepancy between the mass function of structure in molecular 
clouds derived from molecular line observations (shallow), particularly CO (e.g. Kramer et 
al. 1998), and submillimeter-wavelength dust continuum radiation (steep). It is important 
to note that there are exceptions to this trend. Onishi et al. (2002) identified clumps within 
Taurus using H 13 CO + J = 1-0 line emission and found a steep mass function. Also, Kerton 
et al. (2001) and Mookerjea et al. (2004) have found much shallower mass functions associ¬ 
ated with larger, more massive dust continuum clumps in distant Giant Molecular Clouds. 
These latter submillimeter-wavelength structures are most definitely the counterparts to the 
dense cores typically observed in molecular line emission, while the H 13 CO + clumps are 
counterparts to the small-scale, dense dust continuum structures discussed in this paper. It 
remains to be determined, however, whether there is a clear separation between these two 
types of structure or whether the different observing techniques are presently sensitive only 
at the extremes of an underlying smoothly varying distribution. 


4.4. Clump Environment 

The cold dust clumps found are grouped in three main concentrations: B33, NGC2024, 
and NGC2023. These regions were also observed to be bright in line emission; for example, 
in addition to CO and its isotopomers, CS J = 2-1, traditionally a dense gas tracer (LBS). 
Figure 16 overlays the integrated CS observations on the 850 /irn image for comparison. As 
in Orion B North (Paper III), the vast majority of the submillimeter-wavelength clumps he 
above the 3 km s _1 contour, corresponding to column densities of Nh ~ 10 22 cm -2 through 
the CS core region. This is consistent with the result that an extinction threshold may 
exist for the formation of smallscale, thermally dominated, submillimeter clumps (Paper II, 
Paper III, Johnstone, Di Francesco, & Kirk 2005). There are, however, two distinct clumps 
(SMM J054204-02077 and SMM J 054205-02025) residing in particularly low integrated CS 
emission. Both of these sources are highly concentrated but otherwise unremarkable. 

Johnstone et al. (2005) noted that approximately twenty percent of the mass of the core 
region in Ophiuchus has fragmented into identifiable submillimeter clumps. The apparent 
threshold for submillimeter clump formation, however, leads to few sources being detected 
outside the dense core and thus only two percent of the mass of the entire Ophiuchus cloud 
is visible as submillimeter clump fragments. These are intriguing numbers, comparable 
to the fractional mass of embedded young clusters in dense cores (Lada & Lada 2003), 
suggesting that the clumping process may be linked to the formation of stellar associations. 
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Comparing the CS mass determinations of the parent cores, derived by LBS, against the 
total submillimeter clump mass residing within each core for both Orion B South and Orion 
B North (Paper 111) is straightforward. Although there exists a significant spread from 5% 
mass fragmentation in NGC 2071 to 35% mass fragmentation in NGC 2024, overall one 
finds ~ 300 M 0 in submillimeter clumps within ~ 1700 M 0 of CS clumped material, or 
approximately 20 percent fragmentation by mass. This result is identical to that found for 
the Ophiuchus region. 


5. Conclusions 

We have presented images of the southern part of the Orion B region obtained at 
submillimeter wavelengths of 850 and 450 /irn, covering an area 1300 arcmin 2 in extent. 
Utilizing the techniques developed in Papers I, II, and III, we have identified 57 independent 
dust concentrations, or “clumps”, within this region. The majority of these clumps are 
concentrated in two star-forming regions, NGC 2024 and NGC 2023; a third minor grouping 
appears in the Horsehead Nebula cloud. Most of the dust mass is concentrated in NGC 2024. 

We have been able to derive estimated masses and temperatures for the majority of the 
clumps by noting that they can be modelled as Bonnor-Ebert spheres; that is, they can be 
approximated by almost constant density models with low internal velocities. On this basis 
the typical clump temperature is found to be 18±4 K; this value is in reasonable agreement 
with that (21±9 K) derived from the spectral indices obtained from the 850 and 450 /irn 
data. Thus for most of the clumps a temperature of 20 K is a good working approximation. 
However, 9 clumps are centrally concentrated enough that they may be collapsing, and 4 of 
these objects are particularly bright at 850 /mi. All of the latter lie within NGC 2024 and 
have temperatures which are considerably higher than that for the majority of the clumps. 

Using CO isotopomer spectral line data obtained for NGC2023 we have attempted to 
obtain independent estimates of temperatures and masses for the clumps in this particular 
region. With one exception the results obtained agree with those derived from the continuum 
data. It can be difficult, however, to clearly relate CO structures with equivalent continuum 
clumps. This problem becomes more acute at the lower mass end of the clump distribution, 
where continuum sensitivity is limited and individual clumps overlap. 

The mass function of submillimeter clumps found in Orion B South is steep, in agreement 
with the results in other nearby regions (Motte et al. 1998; Paper II; Paper III; Motte et al. 
2001). The fraction of mass in each CS-identified core which has fragmented in observable 
submillimeter clumps is ~ 20 percent. 
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Table 1. Clump properties in Orion B South derived from 850 /jin data. 


Name * 
(SMM J ) 

R . A. b 

( J 2000 ) 

Dec. b 

( J 2000 ) 

Ss50 C 

(Jy) 

crpeakc 

°850 

( Jy / bm ) 

- R e ff C 
( 10 3 AU ) 

C d 

T d d 

( K ) 

log P / k d 
( K / cm 3 ) 

M Td d 

( Mq ) 

Mtj =20 K 

( Me ) 

054086-01528 

05 : 40 : 51.8 

- 01 : 52:48 

1.23 

0.18 

11 

0.40 

18 

6.0 

1.15 

0.98 

054089-02271 

05 : 40 : 53.4 

- 02 : 27:06 

2.18 

0.31 

12 

0.51 

16 

6.0 

2.44 

1.73 

054090-01531 

05 : 40 : 54.0 

- 01 : 53:06 

0.89 

0.21 

9 

0.42 

16 

6.2 

1.00 

0.71 

054091-02278 

05 : 40 : 54.6 

- 02 : 27:51 

2.90 

0.28 

14 

0.45 

19 

6.0 

2.49 

2.31 

054096-02074 

05 : 40 : 57.6 

- 02 : 07:27 

2.18 

0.28 

13 

0.54 

15 

5.9 

2.72 

1.73 

054097-02086 

05 : 40 : 58.2 

- 02 : 08:36 

4.21 

0.40 

16 

0.57 

18 

5.9 

3.92 

3.34 

054100-02092 

05 : 40 : 59.8 

- 02 : 09:15 

5.24 

0.29 

19 

0.52 

19 

5.8 

4.49 

4.16 

054105-01558 

05 : 41 : 03.2 

- 01 : 55:51 

1.43 

0.20 

11 

0.45 

16 

6.0 

1.60 

1.13 

054112-02273 

05 : 41 : 07.0 

- 02 : 27:18 

3.08 

0.48 

13 

0.61 

16 

6.0 

3.46 

2.45 

054116-01580 

05 : 41 : 09.4 

- 01 : 58:00 

1.25 

0.24 

10 

0.48 

15 

6.1 

1.56 

0.99 

054140-01507 

05 : 41 : 23.8 

- 01 : 50:45 

0.96 

0.27 

9 

0.52 

13 

6.1 

1.53 

0.76 

054141-02180 

05 : 41 : 24.8 

- 02 : 18:03 

18.57 

3.19 

21 

0.87 

25 

5.9 

10.72 

14.74 

054141-02289 

05 : 41 : 24.4 

- 02 : 28:54 

0.76 

0.21 

8 

0.44 

15 

6.2 

0.95 

0.60 

054142-02160 

05 : 41 : 25.2 

- 02 : 16:00 

7.42 

0.47 

20 

0.60 

19 

5.8 

6.36 

5.89 

054142-02194 

05 : 41 : 25.4 

- 02 : 19:24 

4.18 

0.48 

14 

0.55 

19 

6.1 

3.58 

3.32 

054146-02208 

05 : 41 : 27.8 

- 02 : 20:51 

2.36 

0.29 

12 

0.45 

19 

6.1 

2.03 

1.88 

054146-02239 

05 : 41 : 27.6 

- 02 : 23:54 

1.90 

0.26 

12 

0.50 

16 

6.0 

2.13 

1.51 

054148-02200 

05 : 41 : 28.6 

- 02 : 20:03 

4.92 

0.48 

15 

0.54 

20 

6.1 

3.91 

3.91 

054149-02212 

05 : 41 : 29.4 

- 02 : 21:15 

5.36 

0.97 

15 

0.77 

17 

5.8 

5.46 

4.26 

054149-02232 

05 : 41 : 29.2 

- 02 : 23:15 

0.83 

0.25 

8 

0.50 

13 

6.2 

1.33 

0.66 

054151-02251 

05 : 41 : 30.4 

- 02 : 25:06 

0.91 

0.19 

9 

0.38 

21 

6.1 

0.68 

0.73 

054152-02262 

05 : 41 : 31.0 

- 02 : 26:12 

3.45 

0.21 

17 

0.45 

19 

5.8 

2.96 

2.74 

054153-02189 

05 : 41 : 31.6 

- 02 : 18:57 

4.56 

0.32 

18 

0.56 

17 

5.8 

4.64 

3.62 

054155-01497 

05 : 41 : 33.2 

- 01 : 49:45 

1.87 

0.32 

10 

0.46 

18 

6.2 

1.74 

1.49 

054156-02209 

05 : 41 : 33.4 

- 02 : 20:54 

1.16 

0.25 

9 

0.47 

15 

6.1 

1.45 

0.92 

054156-02250 

05 : 41 : 33.6 

- 02 : 25:03 

0.78 

0.18 

9 

0.42 

16 

6.1 

0.88 

0.62 

054158-01558 

05 : 41 : 34.8 

- 01 : 55:48 

4.15 

0.40 

15 

0.54 

18 

6.0 

3.87 

3.30 

054158-02303 

05 : 41 : 34.6 

- 02 : 30:18 

0.47 

0.19 

6 

0.39 

16 

6.3 

0.52 

0.37 

054159-01514 

05 : 41 : 35.6 

- 01 : 51:27 

2.01 

0.34 

11 

0.51 

16 

6.1 

2.26 

1.60 

054160-01492 

05 : 41 : 35.8 

- 01 : 49:12 

1.14 

0.41 

8 

0.54 

14 

6.3 

1.60 

0.90 

054160-01564 

05 : 41 : 35.8 

- 01 : 56:24 

15.16 

0.58 

24 

0.53 

26 

5.9 

8.29 

12.03 

054160-02300 

05 : 41 : 36.0 

- 02 : 30:00 

0.49 

0.18 

7 

0.36 

21 

6.3 

0.36 

0.39 

054161-01495 

05 : 41 : 36.6 

- 01 : 49:30 

1.58 

0.44 

9 

0.54 

15 

6.3 

1.97 

1.25 

054162-02172 

05 : 41 : 37.0 

- 02 : 17:15 

19.86 

1.66 

26 

0.82 

24 

5.7 

12.13 

15.76 

054163-01509 

05 : 41 : 37.8 

- 01 : 50:54 

2.18 

0.34 

11 

0.46 

18 

6.2 

2.03 

1.73 

054163-01558 

05 : 41 : 37.6 

- 01 : 55:48 

7.46 

0.58 

15 

0.44 

28 

6.2 

3.69 

5.92 

054164-01504 

05 : 41 : 38.4 

- 01 : 50:24 

1.48 

0.28 

10 

0.44 

18 

6.2 

1.38 

1.17 

054167-02171 

05 : 41 : 40.2 

- 02 : 17:06 

13.91 

1.12 

17 

0.59 

27 

6.2 

7.23 

11.04 

054170-02162 

05 : 41 : 42.0 

- 02 : 16:12 

12.06 

0.68 

23 

0.65 

21 

5.7 

8.91 

9.57 

054172-01543 

05 : 41 : 43.0 

- 01 : 54:21 

85.03 

11.11 

24 

0.86 

47 

6.4 

21.96 

67.50 

054173-01527 

05 : 41 : 44.0 

- 01 : 52:45 

18.04 

0.85 

24 

0.64 

25 

5.8 

10.41 

14.32 

054173-01547 

05 : 41 : 44.0 

- 01 : 54:42 

63.67 

7.99 

23 

0.85 

42 

6.3 

18.82 

50.54 

054173-02170 

05 : 41 : 43.8 

- 02 : 17:03 

7.67 

0.49 

20 

0.60 

20 

5.8 

6.09 

6.09 

054174-01556 

05 : 41 : 44.4 

- 01 : 55:39 

90.41 

15.35 

19 

0.83 

54 

6.7 

19.82 

71.77 

054175-01560 

05 : 41 : 45.0 

- 01 : 56:00 

113.37 

10.15 

29 

0.87 

49 

6.2 

27.86 

90.00 








Table 1 


Continued 


Name a 
(SMM J) 

R.A. b 

(J2000) 

Dec. b 

(J2000) 

Ss50° 

(Jy) 

npeak c 

°850 

(Jy/bm) 

^eff C 

(10 3 AU) 

C d 

T d d 

(K) 

logP//c d 

(K/cm 3 ) 

(M e ) 

Mr d =20 K e 

(Me) 

054182-01576 

05:41:49.4 

-01:57:39 

24.19 

1.07 

27 

0.69 

26 

5.7 

13.23 

19.21 

054182-01596 

05:41:49.2 

-01:59:36 

3.58 

0.47 

14 

0.63 

16 

5.9 

4.02 

2.84 

054185-01534 

05:41:51.2 

-01:53:27 

0.66 

0.19 

8 

0.42 

15 

6.2 

0.83 

0.53 

054188-02002 

05:41:53.0 

-02:00:12 

1.11 

0.30 

8 

0.45 

16 

6.3 

1.24 

0.88 

054190-02004 

05:41:54.0 

-02:00:24 

0.95 

0.31 

8 

0.48 

15 

6.3 

1.18 

0.75 

054192-01565 

05:41:55.0 

-01:56:33 

4.86 

0.22 

19 

0.40 

26 

5.8 

2.66 

3.86 

054195-02009 

05:41:56.8 

-02:00:54 

0.50 

0.25 

6 

0.48 

13 

6.4 

0.80 

0.40 

054204-02077 

05:42:02.4 

-02:07:42 

3.09 

0.95 

12 

0.76 

16 

6.0 

3.47 

2.45 

054205-02025 

05:42:02.8 

-02:02:33 

3.89 

1.14 

11 

0.72 

18 

6.2 

3.62 

3.09 

054207-02017 

05:42:04.2 

-02:01:42 

0.72 

0.19 

8 

0.40 

17 

6.2 

0.73 

0.57 

054207-02032 

05:42:04.0 

-02:03:12 

0.58 

0.21 

7 

0.44 

14 

6.3 

0.81 

0.46 

054218-02042 

05:42:10.8 

-02:04:12 

0.95 

0.30 

9 

0.56 

13 

6.2 

1.51 

0.75 


a Name formed from J2000 positions (hhmm.mmdddmm.m). 
b Position of peak surface brightness within clump (accurate to 3 /; ). 

c Radius, peak flux, and total flux are derived from clfind (Williams et al. 1994). The peak flux and total flux have uncertainties 
of about 20 percent, mostly due to absolute flux calibration. The radius has not been deconvolved from the telescope beam. 

d Quantities derived from Bonnor-Ebert analysis (see text). 

e Mass derived from the total flux assuming T^ = 20 K and «850 — 0.02cm 2 g _1 . 






Table 2. Submillimeter wavelength properties of clumps in Orion B South 


Name 3, 
(SMM J) 

Sitso/Ssso 15 

54 P 50 ak /5 8 P 5 e 0 kb Ti(BE) c 
(K) 

T d { SM) d 
(K) 

054086-01528 

24.4 

23.5 

18 

50 

054089-02271 

8.2 

7.9 

16 

21 

054090-01531 

20.8 

16.6 

16 

50 

054091-02278 

9.4 

8.7 

19 

28 

054096-02074 

10.0 

9.7 

15 

35 

054097-02086 

8.6 

8.3 

18 

22 

054100-02092 

8.6 

8.6 

19 

23 

054105-01558 

12.4 

11.0 

16 

50 

054112-02273 

10.5 

8.7 

16 

43 

054116-01580 

2.9 

5.7 

15 

7 

054140-01507 

8.0 

10.5 

13 

19 

054141-02180 

7.5 

5.7 

25 

17 

054141-02289 

13.1 

12.1 

15 

50 

054142-02160 

8.4 

8.4 

19 

21 

054142-02194 

5.4 

6.5 

19 

12 

054146-02208 

6.9 

7.9 

19 

15 

054146-02239 

11.5 

10.2 

16 

50 

054148-02200 

6.4 

7.2 

20 

14 

054149-02212 

8.0 

6.7 

17 

20 

054149-02232 

8.9 

8.3 

13 

24 

054151-02251 

12.0 

13.0 

21 

50 

054152-02262 

7.9 

9.9 

19 

19 

054153-02189 

2.9 

4.5 

17 

7 

054155-01497 

8.0 

9.2 

18 

19 

054156-02209 

6.1 

7.7 

15 

13 

054156-02250 

13.1 

14.1 

16 

50 

054158-01558 

11.8 

12.5 

18 

50 

054158-02303 

11.6 

11.8 

16 

50 

054159-01514 

13.0 

11.4 

16 

50 

054160-01492 

6.7 

9.9 

14 

15 

054160-01564 

15.5 

19.4 

26 

50 

054160-02300 

13.0 

12.3 

21 

50 

054161-01495 

9.6 

10.4 

15 

29 

054162-02172 

6.8 

5.0 

24 

15 

054163-01509 

13.3 

12.2 

18 

50 






Tabic 2—Continued 


Name 3, 
(SMM J) 

<S'45o/'S'850 b 

Qpeak / Qpeakb 
*^450 /^850 

T d (BE) c 

(K) 

T d (SM) d 

(K) 

054163-01558 

17.5 

18.9 

28 

50 

054164-01504 

10.4 

10.9 

18 

40 

054167-02171 

8.5 

9.1 

27 

22 

054170-02162 

6.0 

8.5 

21 

13 

054172-01543 

14.0 

11.2 

47 

50 

054173-01527 

15.2 

16.5 

25 

50 

054173-01547 

13.9 

11.5 

42 

50 

054173-02170 

6.1 

8.4 

20 

13 

054174-01556 

12.8 

9.8 

54 

50 

054175-01560 

16.3 

12.4 

49 

50 

054182-01576 

22.4 

17.7 

26 

50 

054182-01596 

19.3 

13.9 

16 

50 

054185-01534 

18.5 

24.2 

15 

50 

054188-02002 

17.6 

15.6 

16 

50 

054190-02004 

17.3 

15.8 

15 

50 

054192-01565 

23.4 

29.4 

26 

50 

054195-02009 

15.1 

12.4 

13 

50 

054204-02077 

7.9 

6.9 

16 

19 

054205-02025 

11.3 

9.9 

18 

50 

054207-02017 

13.5 

14.8 

17 

50 

054207-02032 

10.9 

14.3 

14 

50 

054218-02042 

7.6 

8.6 

13 

18 


a Name formed from J2000 positions (hhmm.mmdddmm.m). 

b The 850 /jm and 450 /jm peak flux and total flux are derived from clfind (Williams et al. 
1994) after convolution to an identical beam size. 

c Quantity derived from Bonnor-Ebert analysis (see text). 

d Quantity derived from spectral energy fit to 450 and 850/im integrated fluxes (see text). 






Table 3. 


Dust and gas properties of clumps in NGC2023 


Name 3, 
(SMM J) 

Id (BE) 
(K) 

T d { SM) 
(K) 

Igas(C 18 0) 

(K) 

N(850) b 
(10 22 cm -2 ) 

N(C 18 0) 
(10 22 cm -2 ) 

N(850)/N(C 18 O) 

054141-02180 

25 

19 

18 

12.1 

3.4 

3.6 

054153-02189 

17 

8 

14 

2.1 

4.2 

0.5 

054162-02172 

24 

16 

18 

6.6 

4.9 

1.3 

054167-02171 

27 

24 

25 

3.8 

10.0 

0.4 

054173-02170 

20 

14 

15 

2.6 

5.8 

0.4 


a Name formed from J2000 positions (hhmm.mmdddmm.m). 
“Column density derived assuming the BE temperature for the dust. 
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05 h 42 m 40* 20 s 42 m 00 s 40 s 20 s 41 m 00 s 


R.A. (2000) 


Fig. 1.— The Orion B South region at 850 /mi. North to south, the bright star-forming 
regions are: NGC 2024, NGC 2023, and the Horsehead. Left: 850 /mi emission (greyscale: 
white to black represents -0.25 to 0.5 Jy bm -1 , contours are shown for 0.25, 1, and 4 Jy bm -1 ). 
The named boxes indicate the regions displayed in Figures 4 to 6. The beamsize is noted in 
the lower left of the figure. Right: 850 /mi noise map (greyscale: white to black represents 0 
to 0.15 Jybm -1 ). In the noise map, the lower right hand corner has a slightly higher than 
typical uncertainty of 0.07 Jy pixel -1 while the Horsehead region directly above has a typical 
uncertainty of 0.03 Jy pixel -1 . The greater noise in regions where there is strong signal, such 
as NGC 2024, is due to uncertainties in the calibration of individual bolometers. 
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05 n 42 m 40 s 20 s 42 m 00 s 40 s 20 s 41'"00 s 

R.A. (2000) 


Fig. 2.— The Orion B South region at 450 /urn, covering the same area as described in Figure 

1. Left: 450 /mi emission (greyscale: white to black represents -1 to 2 Jybm -1 , contours at 

2, 8, 32 Jybm -1 ). The beamsize is noted in the lower left of the figure. Right: 450 /.irn noise 
map (greyscale: white to black represents 0 to 1 Jybm -1 ). In the noise map, the lower right 
hand corner reveals a higher than typical uncertainty of 0.9 Jy pixel -1 while the Horsehead 
region directly above has a typical uncertainty of 0.2 Jy pixel -1 . 
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Fig. 3.— Location of the submillimeter clumps found in Orion B South 850 /mi using an 
automated procedure (see text). The circle size approximates the area associated with each 
clump. 
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05 h 41 m 30® 20 s 10' 41"'00' 50' 

RA(2000) 


Fig. 4.— The Horsehead Nebula region (B33). The panels show the 450 and 850 /mi emis¬ 
sion at upper and lower left respectively, plotted on intensity scales chosen to illustrate the 
structure of the emission to good advantage. Overlaid on the 850 //m map are plus-signs de¬ 
noting the locations of the submillimeter sources. At upper right we show the corresponding 
image from the Digitized Sky Survey (DSS), and at lower right an image from the 8-/im data 
obtained from the MSX database. The small circle at lower right in each case indicates the 
respective beamwidths, except for the DSS data. 
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Fig. 5.— The NGC2023 region at submillimeter (left), optical (upper right) and mid-infrared 
(lower right) wavelengths. Other details are as for Figure 4. 
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05 B 4 1 m 56 s 52 s 48 s 44' 40' 36' 


R. A. (2000) 


Fig. 6.— The NGC2024 region at submillimeter (left), optical (upper right) and mid-infrared 
(lower right) wavelengths. The horizontal stripe in the DSS image is due to a bright star 
outside the frame shown here. Vertical banding in the MSX image is an artifact. Additional 
contours for the 450 pm image are 30, 45, 60 and 75 Jy/beam area, and for the 850 pm image 
are 4, 6, 9, 12 and 15 Jy/beam area. Other details are as for Figure 4. 
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Fig. 7.— Derived mass for each of the 57 clumps, assuming Td = 20 K, versus the effective 
radius. Note that the effective radius is derived from the extent of the clump as determined by 
elf ind (see text) and the largest uncertainty in the mass of the clump is the flux calibration, 
which is accurate to about twenty percent. Also plotted (dashed-dotted) are the minimum 
size which a clump might have (resolution limit) and the minimum mass that a clump must 
have for a given radius such that the clump-finding routine can recognize it (approximately 
4 a above the background). The symbols denoting each clump are discussed in Figure 8. 
The three curves from left to right denote the mass-radius relation for Bonnor-Ebert spheres 
with 20 K internal temperatures, an additional equal internal pressure component due to 
turbulence, and external pressures k~ l P = 0.3,1.0, 3.3 x 10 6 Kcm -3 (see text). 
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C 


Fig. 8.— Derived concentration for each of the 57 clumps in the Orion B South region (see 
text). The minimum concentration for a constant density low-mass Bonnor-Ebert sphere 
is C = 0.33, while the maximum concentration beyond which collapse occurs is C = 0.72. 
Clumps with C > 0.72 are denoted by diamonds. (Top) Derived mass of the clump, assuming 
T d 20 K, vs. concentration. (Bottom) Derived radius of the clump vs. concentration. 
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Fig. 9.— Results of determining the physical parameters of clumps using the assumption that 
they are well represented by Bonnor-Ebert spheres with measured concentrations. Plotted 
are the external, confining pressure vs. the internal temperature (assuming an equal contri¬ 
bution from turbulent support). The typical uncertainty in the derived physical parameters 
of an individual clump is similar to the spread in the distribution of points. The symbols 
are the same as in Fig. 8. 




35 



0 10 20 30 40 50 60 

T(BE) (K) 


Fig. 10.— Comparison of the clump temperature derived from the Bonnor-Ebert analysis 
versus the clump temperature derived from spectral energy fitting between the 450 //m and 
850/mi observations (see text). 




13C0 J=2-1 MAP OF NGC 2023 



Right Ascension (J2000) 

Fig. 11.— A 13 CO J = 2-1 raster map of NGC 2023. The beam size is 21". The sampling 
interval is 7.5" in R.A. and Dec. The maximum contour is 58.0 K.km/s and the contour 
interval is 5.0 K.km/s. 
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Cl80 J=2-1 MAP OF NGC 2023 
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Right Ascension (J2000) 

Fig. 12.— A C 18 0 J = 2-1 raster map of NGC 2023. The beam size is 21". The sampling 
interval is 7.5" in R.A. and Dec. The 13 C0 J = 2-1 and C 18 0 J = 2-1 maps were obtained 
simultaneously with a single tuning of the facility A-band receiver. The maximum contour 
is 12.3 K.km/s and the contour interval is 1.0 K.km/s 




Declination (J2000) 



Right Ascension (J2000) 


Fig. 13.— A C 18 0 J = 3-2 raster map of NGC 2023. The beam size is 14". The sampling 
interval is 5" in R.A. and Dec. The maximum contour is 13.4 K.km/s and the contour 
interval is 1.0 K.km/s. 
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R.A.(2000) 

Fig. 14.— The C 18 0 J = 2-1 grayscale image of NGC 2023 (see Figure 12), overlain by 
contours of the 850/jm continuum emision. The contour levels shown are 0.1, 0.2, 0.4, 0.6, 
1.0, 1.5, and 2.0 Jy/beam, and the beamsizes to half-power are indicated at lower right. The 
superposed crosses show the locations of the continuum sources indicated by the automated 
clump-finding routine (see Section 3.1, and also Figures 3 and 5). 
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log(M/M 0 ) 


Fig. 15.— Cumulative number function for the 850 /irn clumps in Orion B South. The thin 
line converts the measured flux directly to mass using a constant temperature Td = 20 K. 
The thick line converts the flux to mass using the derived temperature from fitting the 
clumps to Bonnor-Ebert spheres (see text). The horizontal lines represent the extent to 
which the masses might change due to uncertainties in the Bonnor-Ebert Measurements (see 
Paper III). The steep dashed line has a slope M -1 ' 5 and approximates the high-mass end of 
the Bonnor-Ebert cumulative distribution. The shallow dashed line has a slope M -0,5 and 
approximates the low-mass end of the distribution. The low-mass end is likely to be severely 
incomplete (see text). 
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R.A. (2000) 


Fig. 16.— Contours of CS J = 2-1 emission (from Lada et al. 1991) are overlaid on the 
850 fin i map. The contour levels are at 1.5, 3.0, 6.0, and 12.0 Kkms -1 . Note that a majority 
of the strong dust continuum sources lie above the 3 Kkms -1 contour. 






